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We reproduce the broadband spectrum of Crab pulsar, from UV to very high energy gamma-rays - nearly 
ten decades in energy, within the framework of the cyclotron-self-Compton model. Emission is produced by 
two counter-streaming beams within the outer gaps, at distances above ~ 20 NS radii. The outward moving 
beam produces UV-A-ray photons via Doppler-booster cyclotron emission, and GeV photons by Compton 
scattering the cyclotron photons produced by the inward going beam. The scattering occurs in the deep 
Klein-Nishina regime, whereby the IC component provides a direct measurement of particle distribution 
within the magnetosphere. The required plasma multiplicity is high, ~ 10 6 — 10 7 , but is consistent with the 
average particle flux injected into the pulsar wind nebula. 

The importance of Compton scattering in the Klein-Nishina regime also implies the importance of pair 
production in the outer gaps. We suggest that outer gaps are important sources of pairs in pulsar magneto- 
spheres. 

Cyclotron motion of particles in the pulsar magnetosphere may be excited due to coherent emission of 
radio waves by streaming particles at the anomalous cyclotron resonance. Thus, a whole range of Crab 
non-thermal emission, from coherent radio waves to very high energy 7-rays - nearly eighteen decades in 
energy - may be a manifestation of inter-dependent radiation processes. 

The present model, together with the observational evidence in favor of the IC scattering (Lyutikov et al. 
2012; Lyutikov 2012), demonstrates that the inverse Compton scattering can be the dominant high energy 
emission mechanism in majority of pulsars. 
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1. Introduction 

The pulsar high energy emission is a complicated unsolved problem in high energy astrophysics. It has been been 
under intensive study for nearly four decades (Cheng & Ruderman 1977 Cheng et al. 1986a Romani & Yadigaroglu 



1995 Daugherty fe Harding|1996 Harding et al 



2008 ). The Fermi Gamma-Ray Space Telescope detected a large number 



of pulsars (Abdo 2010b); this revolutionized our picture of the non-thermal emission from pulsars in the gamma-ray 



band from 100 MeV up to about 10 GeV. At even higher energies, in the very-high energy (VHE) band, the detection of 



the Crab pulsar at 25 GeV by the Magic Collaboration ( |Aliu fc MAGIC Collaboration|2008| ) and recently at 120 GeV by 
the VERITAS Collaboration ( |VERITAS Collaboration fe Aliu||2011[ ) in the very-high energy band allow to stringently 
constrain the very-high-energy emission mechanisms 



Geometrical models, based on the idea of the outer gap (Cheng et al. 1986a), are very successful in explaining 
the basic features of the observed 7-ray light curves (e.g., Romani & Yadigaroglu 1995 Harding et al. 2008 Bai & 



Spitkovsky|2010l ). While there seems broad consensus that the particle accelerator is located in the outer magnetosphere, 

Arons||1996 |. 



the radiation physics remain controversial (e. 



Nearly universally, the origin of the emission above ~ 100 MeV was until recently attributed to the curvature 
emission (e.g., Cheng & Ruderman 1977 Cheng et al.||1986a Romani & Yadigaroglu| |1995 Hirotani fc Shibata||1999 
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Hirotani et al. 2003 Hirotani 2007 Harding et al. 20081. For example, Cheng et al. (1986b) concluded that "Crab 



primary outer gap e + /e~ lose most of their energy to curvature 7- rays". Curvature radiation has remained as the 
preferred gamma-ray emission mechanism (Romani 19961 (see also Cheng et al. ( 1986a| 20001; Takata et al. (2008); 

Possible importance of the IC scattering was discussed in application to the Vela pulsar (e.g., 

The IC scattering was 



Tang et al. ( 


2008)). 


Cheng et al. 


1986b; 



1986b| Romani 1996) (see also Cheng et al. 1986a 2000 Harding et al. 2008) 



assumed to be done by the particles in the magnetosphere interacting with surface X-ray photons. The IC scattering 
was not deemed to be the dominant mechanism of high energy emission. 

In contrast, we argued ( Lyutikov et al.|2012 Lyutikov 2012) that the IC scattering may be the dominant source of 
high energy photons in a majority of pulsars (see also Aleksic & et al. 2012). In this paper we further develop the IC 
model to include a modeling of the broadband SED, from UV to very high energy 7-rays, covering nearly ten decades 
in energy. In its essence, the lower energy UV-X-ray peak is due to the cyclotron emission by the secondary particles, 
Doppler boosted by the parallel motion of the plasma to the X-ray range, while the GcV component is due to the 
scattering of the cyclotron photons by the counter-streaming beam, Fig. [I] 



2. Observed spectrum of Crab 



Crab pulsar produces non-thermal radiation from radio to very high energy 7-rays. Radio emission is coherent and 
has a different origin, though it can be related to (or actually trigger) the high energy emission, see Sj6j The non-coherent 
non-thermal emission then spans energies from optical, ~ 1 eV, to very high energy 7-rays, ~ 10 11 eV. We interpret 
the Crab SED as having two spectral bumps (e.g., Kuiper et al.||200T and Fig. [1]), a broad UV-X-ray-soft 7-ray bump, 
~ 1 eV- 10 MeV, and a high energy 7-ray bump, ~ 100 MeV- 100 GeV. In this paper we address the nature of the 
non-coherent non-thermal emission and reproduce the high energy spectrum over nearly 10 decades in energy. 

The low energy part of the Crab pulsar spectrum may be roughly represented as a vF v oc e 1 / 3 for energies below 
~ lOfceV, a flat part between 10 keV- 1 MeV and a falling vF v oc e' 1/3 above MeV (iKuiper et al.||200l|). In the high 



energy band, the spectrum is nearly flat vF v oc e° below few GeVs (Abdo 2010a I and shows a long powerlaw tails 
vF v oc e" 3 8 up to ~ 100 GeV ( |VERITAS Collaboration k. Aliu||201l |). 



Outline of the model 



We assume that emission occurs in the outer gaps with a conical shape region extending between some minimal 
radius r m i n and the light cylinder. The conical shape is naturally an approximation, based on the notion that the high 
energy emission is produced along caustics, where effects of light travel and aberration compensate. We assume that 
the emission region subtends a solid angle Af2. We assume that the gaps have two counter-streaming population of 
secondary particles, each with multiplicity A and constant bulk Doppler factor S. The presence of counter-streaming 
populations of leptons is an important ingredient of the outer gap models, e.g., Cheng et al. (1986a) write "Outer gap 



models will generally result in roughly symmetric flowing streams of relativistic [pairs] in both directions along B field 
lines within the gap". (Halpern & Ruderman 1993 Cheng et al. 1998 also argued that half of the charged particles 



produced gaps move toward the star.) 

In this Section we also assume that particles have constant transverse velocity f3o- Later, in Sj4j in calculating 
the SED this assumption is relaxed. Note that the assumption of constant bulk Doppler factor and constant transverse 
velocity imply violation of the first adiabatic invariant, which would lead to parallel deceleration and increasing transverse 
velocity for incoming beam and parallel acceleration and decreasing transverse velocity for outgoing beam. The implicit 
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Fig. 1. — Overall geometry. It is assumed that the emission region occupies a region of open field lines with solid 
opening angle ~ Af2. The insert shows the radiative processes: two counter streaming beams produce Doppler-boosted 
cyclotron emission and IC emission on the soft photons of the counter-streaming beam. 

assumption here is that there is an external mechanism that leads to the violation of the first adiabatic invariant, see 
Sj6j Thus, we assume that both beams have equal Lorentz factors T and equal densities - theses assumptions can easily 
be relaxed in the follow-up studies. We also assume that we see emission only from the outward propagating beam: the 
direct cyclotron emission and the IC upscattering of the inward propagating cyclotron photons, emitted by the inward 
propagating population. 

Both inward and outward going beams produce cyclotron emission - the cyclotron photons emitted by the outgoing 
beam then produce the UV-X-ray bump. The cyclotron photons emitted by the ingoing beam are IC scattered by the 
outgoing beam producing the GeV emission. Hence, the model can be called cyclotron-self-Compton, CSC below. In 
this Section we approximate the broad observed distributions of UV-X-ray and GeV bumps as having typical values 
e s and tic- Thus, we have four observed quantities (typical energies and fluxes of the cyclotron and IC emission). 
The unknown quantities are the location of the emission zone (local cyclotron frequency) , the relative Lorentz factor of 
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the two plasma components, their density, and the characteristic transverse velocity (which determines the intensity of 
the cyclotron emission). Thus there are four observables and four unknowns. (Additional parameters are the surface 
magnetic field, determined from pulsar spindown, and the beaming solid angle determined from the pulse profile.) Thus, 
neglecting for now the detailed shape of the spectrums, the model has four parameters to be fitted, A,r m i n , /3a, 5. This 
is done using the typical energies and fluxes of the cyclotron and IC bumps . These simplifying assumptions are further 
relaxed in S|4j 



3.1. Kinematics of Compton scattering from counter-streaming beams 

Let us derive the kinematic properties of the cyclotron emission and IC scattering for a system of two counter- 
propagating beams. We assume that the IC scattering is done by the outward propagating particles on the cyclotron 
photons produced by the inward propagating population. As a first step, we assume that there is a typical frequency of 
the cyclotron and IC emission, neglecting the fact that both spectral distributions are broad range. This simplification 
allows us to determine the overall properties of the particles' distribution function required to explain the observations 

Let e denote the photon energy measured in the observer frame (the center of momentum frame) in terms of the 
electron rest mass energy e = e/(m e c 2 ). The cyclotron energy is then measured in terms of the quantum field, Huib ~> b, 
where b(r) = B/Bq, B q = m 2 e c 2 /(eh). 

The forward boosted cyclotron energy is 

e s = 5b^2Tb (1) 

For the backward propagating photon, its energy in the frame of the forward propagating electron is 

€ERF = S 2 b -> 4T 2 b (2) 

where S = + /3)/(l — /?) is the Doppler factor. This will give an IC photon with energy in the lab frame 

c £erf 8r 3 & 
€lC = I+2iE^F TT8f26' (3) 

where the last relations assume F>1. 



The KN transition corresponds to %T 2 KN b ~ 1, 



I 1 ,-1/2 3/2 



TKN = VTb = ^ h ' NS "* (4) 



NS- 



where, b^s = Bns/Bq; for dipolar magnetic field b — 6jvs/%) an d t\r = r/R 

Eqns ([l][3! can be resolved for the local magnetic field at the location of the initial emission of the cyclotron photon 
and the Doppler factor: 

VV e s e/C / [ e s /(2e /c ), e s e IC -> oo 
6 = he (Jl + ^~ ^ JT - K 6IC ° (5) 

where the two limits correspond to Thomson and KN regimes. 
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The rest-frame photon energy is then 



£erf = eic^s [\ 1 + + 1 « i „- - - - ^ ( 6 ) 

In terms of the surface magnetic field and stellar radii, Eqns[5]give for the location of the emission zone 

,1/3 (- ,1/3-1/6 ,-1/2 - - n 

„ Ins _J b NS e ic/ €s ' e » e ic -> u , . 

"•R - . ; ; rr/s ~ 1 r,i/aa/3-i/3,.i/3 . V) 



2/3/ / r~ \ !/ 3 ~ | OVS^/S-I/S/J/S 



The Compton cross-section for forward scattering 

_ 1 + 2M 2 + 2fr 2 <5 4 2 1 + 2l s l IC 2 _ (S-2i IC )(S 2 -2Si IC + 2ij c ) 

" < 1 + 2MJ > 3 rB= ( i+ 2 ilci , { i + ^m^7) )f B= p 

f r| e s e /c -> 



(8) 



where = e 2 /(m e c 2 ) is the classical radius of an electron. The transition to the KN regime occurs at ijc ~ l/(8e s ). 

Let us next apply these relations to the Crab pulsar. For Crab pulsar the surface magnetic field Sjvs = 4 x 10 12 G 
— » bus ~ 0.1. Taking the observed cyclotron peak at ~ 50 keV (e s = 0.1), the IC peak at ~ 1 GeV (ejc = 2000), gives 
the estimates of the minimal emission height and the bulk Doppler factor. 

= 2.5 x 10~ 5 

(9) 



(In Crab the light cylinder is located at r)R t Max 
The KN transition in Crab corresponds to 



VR 


« 15 


6« 


e s 




2e w 


Jw 


4000 


160 


) 



Tkjv « 4 /2 ( 10 ) 



Since 2(5 2 6 sa ^lic^i — 800 3> 1, the scattering occurs in a deep KN regime. 

Both estimates ([9| are reasonable. Most model of the high energy emission place the emission heights at tens stellar 
radii (|Hibschman & Arons| ( |2001[ ); |Baring| ( |2004[ ), Arons, priv. comm.), while radiative model predict bulk Lorentz 



factors in thousands (e.g., 



Hirotani fc Shibata||1999[ ). 



lines, 



3.2. Particle density 

We parametrize the total flux through the gap in terms of the total Goldreich- Julian flux through the open field 

• b N s cR N s 2 

2 r E \c 

(t]q = f2i?jvs/c), times the multiplicity factor within the emission region A, times the relative opening angle of the gap 
region at distance r, Af2c/(7rrfi) = AQ/(jrr]Rr)Q), 

Ar = A^ n ^A^ (12) 

2n T] R r E \ c 
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A number density of particles at a radius r is then 



N, b NS A^ASl cR NS , /1oN 

— « r = -7, ^ — " r (13) 

c 2tt t] R r E X c 



Multiplicity A and location of the emission 77^ are the model parameters to be fitted. 



3.3. Cyclotron fluxes 

In the electron center-of-gyration frame the total single particle photon emissivity (integrated over emission angles, 
photons per second per electron) is 

= 3 ^r^Nfo (e' s ~ bm e c 2 )de' s (14) 
In the observer frame the photon emissivity is 5 times higher 

$h = * I ^f3 2 bS D (e s ~ 5bm e c 2 )de s (15) 
where Ac = h/(m e c) is Compton wavelength, 5]j denotes Dirac delta-function. 



The total cyclotron emissivity is (15) times the number of emitting particles (13): 



1 b 2 Q (3 2 S An A C R NS 
V P h = o i ^3— dr5 D (e - bS) de (16) 

D7T Tj R Ar, 



V P h = —b NS \ m ASll3Z-^de (17) 



Integrating over radius we find the total photon emissivity of a mono-energetic stream of particles: 

1 , . .^^cRj 

Thus a beam propagating in a dipolar magnetic field produces constant luminosity between the minimal and maximal 
limits, e min = b NS 5/if R Max and e max (r) = b NS 5/rf R . 

The minimal energy e m m is emitted at the light cylinder, 

emin = ^ns $ Vn- (18) 

The maximal energy at each location t\r is 

e-min = 5b NS VR^ (19) 
The total cyclotron luminosity, emitted mostly at the innermost limit of the gap, can then be estimated using Eq. 



(|17D as 

Lx = m e c 2 J ~er, ph d~e « ^-b NS X m A^ 2 me< ^ N s i 2 (20) 



Equation (20), together with the estimate of the peak energy, e = Sb^slnn mini gi yes the two constraints on the 
parameters of the model (A, r\R, m im Po, The other two constraints will come from the similar estimates for the IC 
component. 
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3.4. IC emission 

In §3.3| we estimated the cyclotron emission produced by an outgoing flux of particles. To estimate the IC flux we 
assume that the same cyclotron flux is produced by the ingoing beam. Then, at each radius IC scattering occurs on 
the cyclotron photons produced by the ingoing beam between that radius and the maximal radius - the light cylinder 
radius. The total IC emission is then the integral of emissivities over the emission zone, between minimal and maximal 
radii. Thus, for a cold flow, a given eje comes from a cyclotron photon emitted at given height, but the IC scattering 
can occur anywhere inside. 

The maximal IC energy, corresponding to cyclotron emission and IC scattering near the stellar surface, is 

5 3 b NS _ c /0 / 01 n 
e ic = t— ^7 — T2 ~ 8 2 ( 21 ) 
f + 2b NS d' 1 

Where the last relation takes into account that 6jvs ~ 0.1, so that for any 5 > 2 the scattering near the surface occurs 
in the KN regime. 

The minimal IC energy, corresponding to cyclotron emission at the maximal r]R t rnax = cQ/Rns is 



- S3b Ns/V 2 R, mi! 

eW -l + 2b NS Syr,l max 



If the KN regime dominates even for lowest energy cyclotron photons, all the IC photons come out with e/c ~ 8/2 w 7 
independent of the emission height of the target cyclotron photon. 

The photon emissivity rate r\ v h,ic per electron is the number of collisions per unit time of soft photons with the 
electron. The soft photon density in the lab frame, at the location t\r is proportional to the emissivity integrated from 



rj R to the r} Rt Max, Eq. (17 1 



Vph b N S fl2 \ 1 1 (r,o\ 

n ' = ^M^ = -W^ X ^ (23) 

In the forward moving electron frame, the density of soft photons emitted towards the star is T times higher. (Recall 
that we assume that outgoing and inward going particle and photon fluxes are similar). 

n' s = Tn s (24) 

The IC scattering rate in the forward moving plasma frame is then 

Vphje = n 's° c5 ( 25 ) 

In the lab frame it is T times smaller 

Vphjc = Vphjc/ T = n s° c (26) 

Thus, the total IC photon luminosity produced at radius r by the outward moving beam of density Ndr/c is 

dr „ / S 2 e 



Vphjc = n s caN — 8 D [ e IC - - - ^ ) dl IC (27) 



Using photon density (23), number density (13) and integrating over the soft photon energies, using the relation 



h l ~ eic ~ irk) di = 6d ( g - w=ko)) (^fcr (28) 



- 8 - 



the IC photon emissivity (photons per second) becomes 



1 ,2 o2 AO 2 x 2 cr ERj 



Vic = i^birsftAStr&X' ^ a F- eic de IC 



F . _ ^-28i lc + 2ij c \ 



53(6 - 2e/ C ) V^, mm &^ 4 / 3 (<5 - 2e /c ) 2 /3 

see Fig. |4| The maximal IC energy is 



I'tA' iriu. 
■ 2b NS 5' 

The overall shape of the IC spectrum depends on whether a KN regime is reached or not. Far below the KN limit 



1 + 2b NS 5 2 h^ mm 



at energies much smaller than the maxim energy ( 30 ) 



Feic ~ r2 2 (31) 

If the KN regime is not reached, the spectral power smoothly goes to zero as ejc —> tic, max- But if the KN regime is 
reached somewhere in the magnetosphere, the distribution develops a peak near e/c ~ 8/2. The peak is narrower and 
higher in the deeper KN regime: in this case all the photons come out at the energy of the electron, regardless of their 
initial energy. 

The maximal IC energy is 

hc,Max = frjVS-g -^7 To ( 32 ) 



Relations for the soft photon energy and fluxes (^3.3) and for IC luminosity and energy, ( 29p2 ) can be resolved 



for the four parameters of the model A, ?7fl, m m, (3o, 8. Due to the dependence of the cross-section on energy, the general 
solutions of such system are complicated. As we demonstrate below, majority of the IC scatterings occur in the deep 
KN regime. This allows a considerable simplification. 



3.5. Scattering in the deep KN regime 

In the KN limit all the photons come out with the same energy, e/c = 8/2. The cyclotron photon is then emitted 
in a region where b = e s /(2e IC ) -> VR = (2frjvse/c/e~s) 1/3 - 



The IC cross-section in the KN regime is 
The IC emissivity is then 



• - ib (33) 



V P h,ic = ^dr%Po An %^ — 3r8 D {e s - S/2)de IC (34) 

( Z7T As-i 6 S Tlr, 



\ c e s r, R 

Integration over radius gives 



« 1 fh.^X/S 

5/e s )V3 

while integration over soft photon density gives 



: - r ib NS S/e s ) L > \ (35) 



de^ = ^ (36) 

Es £s,Min 
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Fig. 2. — Spectrum of IC photons for mono-energetic beam propagating in the magnetosphere and scattering cyclotron 
photons from a similar inward propagating beam. Different bulk Doppler-Factors are plotted 6 — 100, 400, 4000. For 
higher 6, when scattering enters KN regime, the spectrum shows a peak near e/c ~ 5/2: in this case the photons 
scattered in the KN limit come out at the energy of the electron, regardless of their initial energy. 



The IC luminosity is then 



L 



IC 



-m e c r/phjc 



h 2 

U NS 
1447T 2 



i m e rER\ 



NS j n e s,Ma 



AS 



(37) 



Scaling of the IC luminosity with density squared (oc A 2 ) reflects the fact that both the target photon density and 
the IC scattering rate are proportional to density. The logarithm of the ratio of the highest and lowest energies reflects 



the fact that the soft photon spectral density is constant (see Eq. 17), while in the KN regime the cross-section decreases 
with energy oc 1/e. 



The two expressions for cyclotron (20) and IC luminosity (37) can be resolved for the two remaining parameters of 
the model, the multiplicity A and transverse velocity (3q. Assuming ew5x 10 4 eV/(m e c 2 ) and eje ~ 10 9 eV/(m e c 2 ) and 
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given luminosities Lx ~ 10 35 erg s 1 and Ljc ~ 10 35 erg s 1 , we find 

A = iZL^^^[M — 1 i * 10' 

uns r E m L x 



* = ^fef^ 4x10 ail1 (38) 

We can also verify that the kinetic energy flux, ~ (6 /2)AnGjm e c 3 4:irr 2 Aft is much smaller than the spin-down 
luminosity L$d- This requires 

6 Afl V " ~Y~ = VrAU_ 2 , (39) 



a condition that is well satisfied for parameters in Eq. ( 38 ) . 



Also, we can verify that the ratio of the X-ray luminosity, Eq. ( 20 1 to the the spin-down luminosity, 

_^x_ _ b f lg 3 X6 2 An r ^n- 6 n- 3 - LxX c t e _ 1Q -4 f40) 

L SD ' 32 P ° M A "A C ^ % " b% s m eC 3R% s rt ~ (4U) 

corresponds approximately to the typical observed values (e.g., Kaspi et al.||2006j ). 

The estimates given in this Section demonstrate that under fairly general assumption about the parameters of the 
magnetospheric plasma and, most importantly, assuming a presence of two counter-streaming populations of leptons, it 
is possible to reproduce the overall properties of the pulsar high energy emission. In the following Section we develop a 
more detailed semi-analytical model, taking into account a broad distribution of parallel momenta for both beams. 



4. Broadband model of Crab SED 
4.1. 

In this Section we build a broad-band cyclotron-self-Compton model of the Crab pulsar SED. Note, that the con- 
ventional expressions for the synchrotron-self-Compton (SSC) emissivities that are used, e.g., in studies of blazars (e.g., 



Coppi 1992) are not applicable to pulsar magnetospheres. Typically blazar SSC models assume isotropic distribution 
both of particles and photons in some given frame (blob rest frame moving with a single given velocity) , tangled magnetic 
field and a power law particle distribution function stretching between some minimal and maximal Lorentz factors. All 
of these assumptions are not applicable to pulsars. 

First, in many pulsars (e.g., in Crab), the cyclotron decay times even at the light cylinder are shorter than the 
period. (The spontaneous decay times in the particle rest-frame is 



3 ^ 



(41) 



In Crab, for a particle at rest, the cyclotron decay time is always smaller than a period, r c , by a factor 3 x 10~ 4 at the 
light cylinder. For a relativistically moving particle t c is Doppler-stretchcd; in the observer frame r c = jt' c . It becomes 
of the order of the period at 

BKle^P 1 / 6 300 



VR-c ~ I/!r wrel i/& ~ 7 i/ 6 ( 42 ) 
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Thus, at the light cylinder particles with 7 > 80 are in the slow decaying regime, while at the inner edge of the emission 
region, at rj r ^ m i n — 20, particles with 7 > 10 7 are in the slow decaying regime.) Thus, the distribution function is 
expected to the highly anisotropic, with small pitch angles. Efficient radiative decay may lead to the non-relativistic 
transverse velocities, so that emission occurs not in the synchrotron but in the cyclotron regime. 

Second, the pulsar magnetic field is regular dipolar, changing in its strength over many orders of magnitude within 
the emission region. Also, the cyclotron emission depends on the direction to the observer from a given plasma element. 
Third, and most importantly, the distribution of parallel momenta with pulsar magnetosphere is very broad. Thus, the 
standard off-the-shelf SSC models are not applicable to pulsar magnetospheres. 



4.2. Distribution function 

The distribution function within the pulsar magnetosphere is bound to a complicated, anisotropic function that 
strongly depends on the location within the magnetosphere. A self-consistent model of the high energy radiation 
should take into account evolution of the perpendicular and parallel momenta of particles due to the motion in the 
inhomogeneous magnetic field (e.g., conservation of the first and second adiabatic invariant if applicable), as well as 
their evolution due to the other EM processes that can excite or de-excite the transgression motion and/or the parallel 



component of the momentum (e.g., excitant of transverse motion due to the anomalous Doppler resonance, see Kazbegi 



et al. (1991); Lyutikov et al. (1999a) and fj6] It is currently beyond our abilities to calculate its details and its evolution 



within the magnetosphere reliably and self-consistently. We can only hope to catch its main properties by appealing 
to the basic theoretical ideas and using observations to probe it. The theoretical assumption will naturally give only 
an approximation, hopefully self-consistent within a given model. The requirement on the model then is that using the 
minimal number of parameter it should be able to reproduce the overall broadband properties of the pulsar emission. 
Below we describe such a model. We stress that due to the complicated nature of the problem, the fit parameters we 
derived for the particle distribution are probably not unique and also, and not precisely determined within this simple 
model. Yet, the model is able to reproduce the bulk properties of the Crab high energy emission using one major 
assumption (that of a presence of counter-streaming plasma components) and a number of parameters derived from 
observations, like the particle spectrum and its evolution with radius. 

As mentioned above, we assume that the gaps have two counter-streaming population of secondary particles. 



This has been a common assumption in many models, e.g., Cheng et al. (1986a). We also assume that we see emission 



only from the outward propagating beam: the direct cyclotron emission and the IC upscattering of the inward propagating 
cyclotron photons, emitted by the inward propagating population. This produces two spectral bumps, in the UV-X-ray 
range, which we would call the low energy bump, and the > 100 MeV - GeV feature, which we would call the high 
energy bump. 



4-2.1. Parallel distribution 



Different parts of the spectrum probe different parts of the particle distribution. Since in the KN regime of the 
IC scattering the energy of the photons is of the same order as the energy of the scattering electron regardless of the 
target photon energy, the high energy bump provides a direct measurement of the bulk particle population (weighted 
by the energy dependence of the KN cross-section). Since above break the observed spectral index is < —4 (VERITAS 
Collaboration fc Aliu||2011 ), this required that the particle spectrum above the break is < — 3 (in the deep KN regime). 



Below the break the the observed spectral index is flat Abdo (2010a), this required that the particle spectrum well below 
the break is ~ 1. The modeling also then requires that there should be another energy region with the particle spectrum 
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below the break is ~ —1. 

The high energy part of the low energy bump provides another constraint on the distribution function at the 
highest energies. Note that if at the highest energies the particle spectrum is a power-law, / oc 8 , the spectral energy 
distribution is oc e 2 i] s oc e 3 ~ l . Thus, to have a finite emissivity one either needs I > 4 or there should be an upper cut-off 
to the particle distribution. (For \l\ < 3 the cyclotron component keeps rising at high energies, contrary to observations). 
Also, there should be an exponential cut-off to the particle energies to avoid cyclotron component showing up in the 
VERITAS band. 

In addition, there should be a minimal Lorentz factor 70, of the order 7 m i„ <~ l/\/Af2 ~ 10. In this case the 
observed profile is mostly determined by the geometrical factors and is nearly energy-independent. (In addition, our 
simplifying assumption of the IC scattering in the deep KN regime is likely to break down for the low energy tails of 
the distribution function, see Eq. Q.) 

Summarizing, the following distribution function of parallel momenta is inferred from observations, see Fig. [3j 



/(*)« 



5 m , m « 1, for S Q < 5 < Si 

S n , n rj -1, for Si < S < S 2 

S l , I w -3, for S 2 < S < S 3 
exp S/S Sj for8 3 <S 

S ~ 10, Si re 10 4 , S 2 « 10 6 , S 3 > 10 8 (43) 

(The value of the exponential cut-off is not well determined, 10 8 < ^3 < 10 12 ). 

As the estimate of the highest energy that a particle can reach, we can balance curvature losses in the magnetic 
field with the curvature radius of the order of the light cylinder, R c = QRlc, with acceleration by the electric field of 



the order of the magnetic field, Lyutikov et al. (2012) 



Note, that the derived 73 « 2.5 x 10 8 (see below) is close to this theoretical limit (for smaller 83 the cyclotron spectrum 
has an upper cut-off below ~ MeV energies.) 

Most importantly, we assume that the parallel distribution function is constant throughout the magnetosphere - 
an obvious simplifications. This also implies that the motion of particles in the magnetosphere is non-adiabatic. For 
particles gyrating in the dipolar magnetosphere, and in the absence of any other interaction, there an effective parallel 
and perpendicular force due to the conservation of the first adiabatic invariant. The assumption of non-adiabaticity 
(besides being a highly simplifying in the model) implies (and is justified by the assumption) that other non-adiabatic 
forces are at play: particles are accelerated by the parallel electric fields, they produce secondary pairs via various 
radiative processes, are excited to higher Landau level by absorbing radiation, simultaneously lose their transverse 
energy by emitting cyclotron photons and, possibly gain transverse energy by emitting photons at the anomalous 
Doppler resonance, see Lyutikov et al. (1999a I and fj6] 



In addition, we neglected the influence of the curvature and IC radiation reaction on the particle motion. To verify 
the validity of this assumption let us compare a total electron flux N e , normalized to the GJ flux through open field 



1 In this context we note, that in principle all the pulsar high energy emission, from X-ray to 7-rays, can be modeled as boosted cyclotron 
emission with a population of particles with I ~ 3. The GeV bump is then just a bump in the parallel distribution. Variations in the X-ray 
and 7-ray profiles then can be due to somewhat different spacial distributions of corresponding particles. We disfavor this possibility since 
such a model would require highly correlated X-ray to 7-ray signals. 
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Fig. 3. — The parallel distribution function f(S), Eq. (43). In addition, there is a counter-streaming beam with the 
same distribution. Various parts of the SED constrain various parts of the distribution function; see text for details. 
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(45) 



Thus, for A >> 10 6 each secondary particle emits more than one photon - radiative drag is then important in changing 
the parallel momentum (since in the KN regime each scatter results in a particle's momentum change of the order of 
unity). Since our model predicts A ~ 10 6 (see below), as a simplifying assumption we can neglect radiation reaction on 
the parallel particle motion. 



4-2.2. Transverse distribution 



Within the CSC model there are two unknown functions: distributions of parallel and perpendicular momenta. For 
the perpendicular distribution we assume that in the center of gyration of the emitting leptons the transverse velocity 
is non-relativistic, so in that frame the photons are emitted at the local cyclotron frequency. Then, the transverse 
distribution function controls mostly the intensity of the cyclotron emission. In principle the value of the transverse 
velocity should be calculated self-consistently using some excitation mechanism, balanced by radiative losses. A possible 
excitation mechanism is outlined in f|6] For now, we just use a parameterization fio(r). We note that the transverse 
velocity should be negligible below some radius r m i n (otherwise the cyclotron photon peak frequency will be at too high 
frequencies). Somewhat arbitrarily we chose /3o(r) — fio(rj r — Tj^min) 3 (see, though, ^6]for a possible justification of this 
scaling). Partly this choice is motivated by simplicity and the post-factum good resulting fit. Also, this choice of the 
perpendicular velocity is not independent, it is a function of the parallel velocity distribution. 



4.3. Cyclotron emission 

In f|3]we presented simple order-of-magnitude estimates that demonstrate the validity if the CSC model. Let us next 
add a spread in parallel momenta and reproduce some of the details of the observed spectrum. The above relations for 
cyclotron emission produced by the mono-energetic beam can be easily generalized to a broad distribution by multiplying 



the practice density (Eq. (16 1 by f(S)dS , N — > Nf(S)dS and integrating over the Doppler factor S. (For convenience 



we express the distribution function in terms of the Doppler factor 8 and not the particle momentum.) 
The cyclotron emissivity for a broad distribution is then 

Vs = ^-b 2 NS p 2 Q 6An m \ Cl ^6 D (e - b NS 6/ri 3 R ) f(S)dS ^ de (46) 

A c 7 1r 

Integrating over 8 we find 

r,, = ^ 2 An m X C -^f (g) ~ed~er, R dr, R (47) 



To proceed further (the integration over radius) we assume a particular distribution function given by Eq. (43) The 



photon spectral density n s (e,rj R ) at each radius rj R is then found by integrating the cyclotron emissivity (47) from this 
radius to the maximal (light cylinder) radius. The total cyclotron flux is then found by integrating the emissivity (47) 
from some minimal radius (parameter r\ r , m in to be fitted) up to the light cylinder (parameter rj r , max = 160). 



4.4. IC emission 

After the soft photon spectral density is found at each radius we calculate the IC emission. Calculations of the IC 
emission can be done in a general case, but since the scattering is generally in the deep KN limit, the corresponding 
relations simply considerably. In this case the IC photons come out with distribution oc 8d{^ic ~ 8 /2)f($)d8, where 
f(8) is the distribution function of parallel momenta. If the density of the soft photons at radius r/ R is n s {e s ,rj R )dt Sl 
the IC luminosity is 

mc = ^Ml m \n a (i.,T, R )^8 D (6 - 2~e IC ) C ^^f(S)dSde IC (48) 
47r e s Ac r\ R 
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Integration over 5 gives 

e s VR 

A = ^Afi^A 07 ^ (49) 

Coefficient A is the overall normalization, while the shape of the IC spectrum is determined by the particle distribution 
function /, convolved with the photon density n s , dependence of the KN cross-section, oc l/e s and the dependence of 
the soft photon density on radius. 



4.5. Fit to data 



The model is compared with the data in Fig. Q . We stress that this is only a comparison and not a multi-parameter 
fit. Thus, we do not plot the error bars, nor calculate various goodness-of-fit quantities. The purpose of this figure 
is not to make a detailed fit of the spectrum (this would require a more detailed modeling), but to demonstrate that 
under very general assumption the CSC model reproduces main observational features. For this particular figure the 
parameters are Sq = 1.5 (the minimal energy cut-off; the value of Sq determines the low rise energy of the IC component), 
m = 1 (particle energy spectrum below the first break; m is determined by the IC spectrum below the break, in the 
Fermi band, (Abdo 2010a)), Si — 3 x 10 4 (the first particle break energy), n = — 1 (particle energy spectrum above 



the first break), 62 = 5 x 10 6 (the second particle break energy; both n and 62 determine the higher energy part of the 
cyclotron SED), I = —2.9 (particle energy spectrum above the second break; I is most sensitive to the IC spectrum above 
the break determined by VERITAS Collaboration & Aliu (2011)), S 3 = 5 x 10 s (the exponential cut-off), rj r ^ min = 20 



(the minimal distance from the star where cyclotron emission is produced, A = 5.5 x 10 6 (average particle multiplicity; 
A determines the relative scaling of the cyclotron and IC component), (3 = 9 x 10~ 1Q (r] r — r) r _ min ) 3 (the transverse 
velocity in the center of momentum frame; /3q determines the overall normalization, scaling with radius, chosen by the 
prediction of the balance between quasilinear diffusion at the anomalous resonance and spontaneous emission at the 
normal cyclotron resonance, see mostly determines the shape and the relative intensities of the low and high energy 
tails of the cyclotron component). (The small value of So formally violate our assumption of the IC scattering in the deep 
KN regime, but we deem it acceptable giving the simplicity of the model and the uncertainties of the multi-parameter 
fit.) 

The required plasma multiplicity, A ~ 10 6 — 10 7 , is higher than is typically achieved in pair production models 
(Arons & Scharlemann 1979 Hibschman & Arons 2001 Harding fc Muslimov||20flj , but is comparable to the average 
multiplicity of A ~ 10 6 required by the observations of the Crab nebula's X-ray emission (Shklovsky 1970; Arons|l9"96 l, 
see also Sj5] 



5. Two photon pair production in outer gaps 



The models of the pulsar magnetospheres are mostly based on the Ruderman & Sutherland ( 1975 1 model of pair 
production near the polar caps. It was then expected that polar cap regions are intense sources of high energy emission 
( paugherty fc Harding| 19821 11996| |Harding fc Muslimov|1998{ |Zhang & Harding|2000| |Baring|2004) . These expectations 
are not supported by few years of Fermi data, no gamma-ray emission from pulsar polar caps have been observed. This 
calls into question the very paradigm of polar cap particle creation. Previously, pair production in the outer gaps 



between the 7-ray photons and surface X-ray photons were considered by Cheng et al. ( 2000 ) 
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Fig. 4. — The broadband spectrum of the Crab approximated with the CSC model. The data are from |Kuiper et &L\ 



(2001); Abdo (2010a); VERITAS Collaboration & Aliu (2011). The IC bump in the KN regime provides a direct 



measurement of the bulk particle distribution, while the high energy part of cyclotron bump constrains the very high 
energy tail of the particle distribution. This is a fit over nearly ten decades in energy, using only a handful of parameters. 



On the other hand, the Crab nebula X-ray emission requires a huge flux of particles, with the average over the 
open field lines multiplicity of < A >^ 10 6 (Shklovsky 1970). Thus, there is a contradiction between polar cap pair 
production models and observations. The current model provides an interesting new possibility for pair production. IC 
scattering in the KN regime is tightly related to the two photon pair production: the corresponding cross-sections are 
very similar (e.g., |Blumenthal fc Gould||1970[ [Anaronian||2004| |Dermer fc Menon||2009[ ). 

The two photon pair production on photons of energy e\ and 1i is a threshold process, which requires sq — e\t2 > 1< 
For a soft photon produced at region of magnetic field b (in terms of quantum magnetic field) and having reduced 
energy in the observer frame <L\ = 5b, interacting head-on with the IC photon scattered in the KN regime, ea ~ 6/2, the 
threshold condition reads 



65 2 /2 = l^ = (W2) 1/3 <5 2/3 



100 



(50) 
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for i5 ~ 2ejc ~ 4000. Thus, the threshold conditions for two photon pair production are satisfied nearly in the whole 
magnetosphere of Crab pulsar (recall that in Crab the light cylinder is located at r)R = 160 and both the soft and the 
VHE photons have a broad spectral distribution). 



The pair-production cross-section has a maximum at the level of <r 77 « Q.2ctt achieved at sq ~ 3.5 — 4 ( Aharonian 
2004] ). For s > 1, 

<r 77 « (2/3)<tt 1 ^ S0 (51) 



so 

„2~ ™ „2\ 



Estimating the photon density as n p h ~ Lx /(Aflr Ijcm^c \ the optical depth for the two photon pair production at 



the location ( 50 lis very high: 



r 77 « m 77 n ph « -^3 -j^- — — = 10 Afi 2 (52) 

Such a very high optical depth is, naturally, an overestimation, since then no 7-ray emission would then be expected 
(c/. the model of gap closure by the photon-photon pair-creation process between the high-energy 7-rays emitted in the 
gap and the X-rays coming from the stellar surface by Zhang & Cheng 19971. For example, in the current estimate we 



assumed that all the soft and VHE photons interact in head-on collisions at a peak of pair production cross-section - 
clearly an upper limit. Still, a very large value of t 77 implies that a 7 — 7 pair production between the soft and VHE 
photons is an important process in pulsar magnetospheres. 

The possibility that outer gaps are the main sources of pair leads to the requirement of the very high particle 
density. If the average over the open field lines multiplicity is < A >^ 10 6 , since the gaps occupy a small fraction of 
the open field lines solid angle, the multiplicity within the gaps should be A ^< A > tttj^/AH, w 2 < A > Afl^VB,', i-c, 
tens to hundreds times higher. This compares favorably with our independent estimates, §4.5| 



To construct a self-consistent model of pair production in the outer gaps is a formidable non-linear problem: the 
plasma density depends on the photon field, which in turn depends on the plasma density; the accelerating electric field 
also depends on the plasma density; the excitation of particle gyration, which controls the cyclotron emission rates, 
also depends on the details of the distribution function (both parallel and transverse, via the growth rate of the maser 
instability), while the evolution of the distribution function also depends on the state of particle gyration, see next 
Section. 



6. Relating radio to high energy emission 



The present model of pulsar high energy emission requires that the emitting particles have a finite pitch angle. This 
requires a mechanism that would excite particle gyration. In this Section we outline such a model, to be addressed in 
more details in a subsequent publication. 

A possible excitation mechanism is related to the generation of the pulsar radio emission at the anomalous cyclotron 
resonance ( Kawamura fc Suzuki|1977 Machabeli fc Usov|1979 Kazbegi et al.|1991 Lyutikov et al.|1999a ) . In this model 
the pulsar radio emission is produced directly by maser-type plasma instabilities operating at the anomalous cyclotron- 
Cherenkov resonance uj — kuvu + |wg|/ Ires — (note the sign in front of the cyclotron term). The instabilities are due 
to the interaction of the fast particles from the primary beam and the tail of the distribution with the normal modes of a 
strongly magnetized one-dimensional electron-positron plasma. The waves emitted at these resonances are vacuum-like, 
electromagnetic waves that may leave the magnetosphere directly ( Lyutikov et al.|fl999b ). 

At the anomalous resonance a particle emitting a wave undergoes a transition up in Landau levels. Thus, initially 
one-dimensional distribution develops a finite pitch angle. As a result, resonant particles start emitting cyclotron 
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photons at the normal cyclotron resonance. The transverse particle distribution can achieve a balance between the 
diffusive spread due to the coherent interaction with the waves at the anomalous resonance and spontaneous photon 
emission at the normal resonance ( Kazbegi et al.|[T991| Chkheidze et al.||2011 Lyutikov||1998b Chkheidze et al.||2011 |. 
In our notation, the typical transverse momentum (Eq. (11) of Chkheidze et al.||2011 ) is 

a ~ 1 ^ c 3 3/4 

dns far^s) ' 

The scaling with radius is the one chosen for the fit, §4.5| Typical values of /3q derived using quasilinear diffusion, Eq. 



(53), are of the same order of magnitude as required by the spectral fit, £4.5 



In addition, other radiative processes will affect the particle distribution function: (i) excitation of cyclotron motion 
due to IC scattering and (ii) during two photon pair production; (iii) induced Raman scattering within the magneto- 
sphere, ( Lyutikov||1998a l; (iv) single particle cyclotron absorption (Petrova & Lyubarskii 2000). Since these processes 
can be nonlocal within the magnetosphere (a photon emitted at one location can be absorbed/scattered at a very 
different location), this presents a complicated electrodynamic problem. 

In this Section we discussed how excitation of coherent radio wave at the anomalous cyclotron resonance in the 
outer gaps can excite the cyclotron motion of the resonant particles. Spontaneous emission of the cyclotron photons 
then produces the UV- X-ray bump, while IC scattering produces the VHE emission. Due to the highly relativistic 
motion all the emission components are beamed along the local magnetic field: this explains the similar profiles of Crab 
pulsar from radio to VHE gamma-ray emission. Thus, the model unites the Crab pulsar non-thermal emission from 
tens megahertz to hundreds of GeVs, nearly eighteen decades in energy. Note that the outlined model, uniting radio 
and gamma ray emission, is applicable specifically to Crab pulsar, and specifically to the main pulse and the interpulse, 
and not the radio precursor. Recall, that Crab is one of the few pulsars where radio and high energy peaks are nearly 



aligned in phase. In most pulsars, radio peak precedes the gamma rays peaks (Abdo 2010b). A relatively weak radio 
precursor in Crab is commonly associated with the "normal" radio emission, while the main pulse and the interpulse 
are, probably, "different" emission mechanisms {e.g., Moffett k. Hankins|[l~996 ) . 



7. Conclusion 

In this paper we demonstrate that the cyclotron-self-Compton model is a viable model of pulsar high energy emission. 
Using observationally-constrained properties of the distribution function of particles within the pulsar magnetosphere, 
as well as fairly general additional theoretical assumptions, we are able to reproduce the overall SED of the Crab pulsar 
over nearly ten decades of energy. The key theoretical assumption is the presence of two counter-streaming populations 
within the outer gaps. The outward moving beam then produces UV- X-ray photons via Doppler-booster cyclotron 
emission, and the GeV emission by Compton scattering the cyclotron photons produced by the inward going beam. As 
a simplifying assumption, the parameters of the inward and outward going beams were chosen to be the same: if the 
distribution is skewed in favor of the outgoing beam, the intensity of the ingoing radiation will be reduced, producing 
a nearly one-sided emission pattern (in addition, cyclotron absorption may be important for the inward going radiation 
( |Ruderman|[T995| |Cheng et al.p000| |). 



Naturally, the current one-dimensional model is a simplification; it is expected that more advanced models, that 
incorporate self-consistently the geometry, structure of accelerating electric fields and various radiative processes, will 
result in the modification of the inferred parameters. The most important simplifications include, first, the assumption of 
a IC scattering in a deep KN regime is not valid for scattering of the cyclotron photon emitted close to the light cylinder 
by the slowly moving particles. Second, the separation of the distribution function into parallel and perpendicular parts 
is likely not to be a unique from the observation point of view. For example, relativistic transverse motion vill increase 
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the emitted frequency in the center-of-gyration frame, thus reducing the required parallel boost. 

I would like to thank Jonathan Arons, Charles Dermer, Kouichi Hirotani and George Machabcli for encouraging 
discussions. 
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